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ABSTRACT 


We derive the primordial power spectra and spectral indexes of the density 
fluctuations and gravitational waves in the framework of loop quantum cosmol¬ 
ogy (LQC) with holonomy and inverse-volume corrections, by using the uniform 
asymptotic approximation method to its third-order, at which the upper error 
bounds are < 0.15%, and accurate enough for the current and forthcoming cos¬ 
mological observations. Then, using the Planck, BAO and SN data we obtain 
the tightest constraints on quantum gravitational effects from LQC corrections, 
and hnd that such effects could be well within the detection of the current and 
forthcoming cosmological observations. 

Subject headings: Loop quantum cosmology, inflationary power spectra. Uniform 
asymptotic approximation, precision cosmology 
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Introduction 


Quantization of gravity has been one of the main driving forces in physics in 
the past decades ([Kiefer 2012), and various approaches have been pursued, including 


string/M-Theory (Becker et ah 2007), loop quantum gravity (Rovelli & Vidotto 2015), and 


more recently the Horava-Lifshitz theory (Hofava 2009). However, it is fair to say that our 


understanding of it is still highly limited, and none of the aforementioned approaches is 
complete. One of the main reasons is the lack of evidences of quantum gravitational effects, 
due to the extreme weakness of gravitational fields. 

This situation has been dramatically changed recently, however, with the arrival of 


the era of precision cosmology (Kiefer & Kramer 2012 Krauss & Wilczek 2014; Woodard 


2014). In particular, cosmic inflation (Guth 1981), which is assumed to have taken place 
during the first moments of time, provides the simplest and most elegant mechanism 
to produce the primordial density perturbations and gravitational waves. The former 
is responsible for the formations of the cosmic microwave background (CMB) and the 


large-scale structure of the universe (Baumann 2009). Current measurements of CMB 


Komatsu et al. 

201l| 

BICEP2/Keck and Planck Collaborations 

2015 

Planck Collaboration 


2015) and observations of the large scale distributions of dark matter and galaxies in the 


universe (Eisenstein et al. 2005 Tegmark et al. 2006; Beutler et ah 2011 Blake 2011) are 


in stunning agreement with it. On the other hand, since inflation is extremely sensitive 


to the Planckian physics (Baumann 2009 Burgess et al. 2013 Brandenberger & Martin 


2013), it also provides opportunities to get deep insight into the physics at the energy 


scales that cannot be reached by any of man-made terrestrial experiments in the near 
future. In particular, it provides a unique window to explore quantum gravitational effects 
from different theories of quantum gravity, whereby one can falsify some of these theories 


^The corresponding author 
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with observational data that have the uncomprehended accuracy (Abazajian et ah 2015), 


and obtain experimental evidences and valuable guidelines for the final construction of the 
theory of quantum gravity. 

In this Letter, we shall study the quantum gravitational effects of LQC in inflation 


(Bojowald 2005 Ashtekar & Singh 2011 Barrau et al 2014), and show explicitly that these 


effects could be well within the detection of the current and forthcoming cosmological 


experiments (Abazajian et al. 2015). Such effects can be studied by introducing appropriate 


modihcations at the level of the classical Hamiltonian, very much similar to those studied 


in solid state physics (Bojowald 2005 Ashtekar & Singh 2011; Barrau et al 2014). It was 


found that there are mainly two kinds of quantum corrections: the holonomy (Mielczarek 


2008a, 2009 Mielczarek et al. 2010 Grain et al. 2010 Li fc Zhu [2011 ; Mielczarek et al. 


2012 Cailleteau et al. 2012b|[a Mielczarek 2014), and inverse-volume corrections (Bojowald 
&: Hossaln 2007 2008a|[b ; Bojowald et al. 2009 2010 Amoros et al. 2014 ; Bojowald &: 
Calcagni |2011 ; Bojowald et al. 2011a|[b ). These corrections modify not only the linear 
perturbations, but also the space-time background. 

In particular, for a scalar field 0 with its potential V{(j)), the holonomy corrections 
modify the Friedmann and Klein-Gordon equations to the forms. 


^2 ^ SirGa^p^ 


^ P<f> 


^ \ pc 

0" + 2-H0' + = 0, 


( 1 ) 

( 2 ) 


where a denotes the expansion factor, "H = a'/a, and a prime denotes the derivative 
with respect to the conformal time rj (= J dt/a{t)). pc is a constant and characterizes 
the energy scale of the holonomy corrections, with p^ = 0'^/(2a^) -|- V{(j)). Clearly, the 
big bang singularity normally appearing at p^p = oo now is replaced by a big bounce 
occurring at pp = pc- In the infrared (IR) we have ppjPc 1, and Eq.([^ reduces to that of 
general relativity (GR). The evolutions of the anomaly-free cosmological scalar and tensor 
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perturbations are described by the mode function /ifc(? 7 ), satisfying the equation (Cailleteau 
et al. 2012a[b ), 


f^kiv) + ( ) hfc(h) = 0, 


( 3 ) 


where ojUrj) = Vt{rj)k‘^ with klirf) = 1 — ‘Ipj Pc- The background-dependent function z{p) is 
given by zs (= acp'/'H) for the scalar perturbations, and zt (= a/y/U) for the tensor ones. 
To the first-order of the slow-roll parameters and 6h{= p/Pc "C 1), the inflationary spectra 
and spectral indexes with the holonomy corrections have been recently obtained, by further 


assuming that the slow-roll parameters and 5h are all constants (Mielczarek 2014). 


With the inverse-volume corrections, on the other hand, the Friedmann and 
Klein-Gordon equations are modified to the forms (Bojowald & Calcagni 2011[), 




( 4 ) 

( 5 ) 


in which p = a?, a = 1 + aoSpi + O (Spi), "d = 1 -|- ’do6pi -|- O {Spi), and 6pi oc a where 


ao, 'do and a are constants [Note that here we use'd instead of z/ adopted in (Bojowald 


& Calcagni 2011), and reserve u for other uses.]. The values of ao and a are currently 
subject to quantization ambiguities, while the magnitude of 6pi is unknown, as so far we 


have no control over the details of the underlying full theory of quantum gravity (Bojowald 


& Calcagni 2011). However, when a takes values in the range 0 < a < 6, the size of 6pi 


does not depend on ao and "do, and can be written in the form 6pi = {api/aY, where api 
is another arbitrary constant. The constant 'do is related to Oq and a via the consistency 
relation 'do(cr — 3 )(ct -|- 6) — 3ao(cr — 6) = 0. However, to make the effective theory viable, we 
shall assume Spi{p) -C 1 at any given moment, so we can safely drop off all the second- and 
high-order terms of 6pi{p). This assumption also guarantees that the slow-roll conditions 
can be imposed, even after the inverse-volume corrections are taken into account. 
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With the above assumption, Bojowald and Calcagni (BC) (Bojowald & Calcagni 2011) 


studied the scalar and tensor perturbations with the inverse-volume corrections, and found 
that the corresponding mode function Hkiv) can be also cast in the form ([^, but now with 

ujUt]) = {l + 2ao6pi{7]))k‘^ ( 6 ) 


for tensor, and 


^kiv) = (I + 2f3o5pi{ri)) 


(7) 


for scalar, where (3o = cri9o (a -|- 6 ) /36 -|- ag (15 — a) jVl. With such modihed dispersion 
relations, BC calculated the corresponding power spectra and spectral indexes to the 
hrst-order of the slow-roll parameters, from which, together with Tsujikawa, they found 
(Bojowald et ah 2011a||b ) that the LQC effects are distinguishable from these of the 
noncommutative geometry or string, as the latter manifest themselves in small scales 


(Tsujikawa et ah 2003 Piao et ah 2004; Calcagni & Tsujikawa 2004), while the former 


mainly at large scales. To hnd explicitly the observational bounds on the inverse-volume 
quantum corrections, they considered the CMB likelihood for the potentials V(</>) = Xn4>^ 
and V (0) = Voe~^^'^, by using the data of WMAP 7yr together with the large-scale 
structure, the Hubble constant measurement from the Hubble Space Telescope, supernovae 


type la, and big bang nucleosynthesis (Komatsu et ah 2011 Reid et al. 2010| Riess et al 


2009 Kowalski et al. 2008 Buries & Tytler 1998), the most accurate data available to 


them by then, and obtained various constraints on 6{k) for different values of a at the 
pivots kg = 0.002 Mpc“^ and kg = 0.05 Mpc“^ where 6{k) = ag5pi{k) for cr 7 ^ 3, and 
5{k) = dg5pi{k) for cr = 3. An interesting feature is that the constraints are very sensitive 
to the choice of the pivots kg, specially when a is large (cr > 2 ), but insensitive to the forms 
of the potential V{(j)). 


In this Letter our goals are two-fold: First, we calculate the scalar and tensor power 
spectra, spectral indexes and the ratio r to the second-order of the slow-roll parameters, for 
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both of the holonomy and inverse-volume corrections, so that they are accurate enough to 


match with the accuracy required by the current and forthcoming experiments (Abazajian 


et ah 2015). This becomes possible, due to the recent development of the powerful uniform 


asymptotical approximation method ( |Habib et ah 2002; Zhu et ah 2014a|[b|[cp! ), which 
is designed specially for the studies of inflationary models after quantum gravitational 
effects are taken into account. Up to the third-order approximations in terms of the free 
parameter (A“^) introduced in the method, which is independent of the slow-roll inflationary 


parameters mentioned above, the upper error bounds are less than 0.15% (Zhu et al 


2014d). Second, we shall use the most recent observational data to obtain new constraints 


on 6{ko) for the power-law potential U(0) = where n is chosen so that r < 0.1. 

With such constraints, we shall prove explicitly that the quantum gravitational effects from 
the inverse-volume corrections are within the range of the detection of the forthcoming 


experiments, specially of the Stage IV ones (Abazajian et al. 2015). 


2. Inflationary Spectra and Spectral Indexes 


To apply the uniform asymptotic approximation method, we first rewrite Eq.([^ 
to = [X^g{y) + (l{y)] i^k{y) I where y = —kr], and the parameter A is a large 

constant to be used to trace the order of approximations. The reason to introduce two 
functions g{y) and q{y), instead of only one, is to use the extra degree of freedom to 
minimize the errors ( [Zhu et ah|2014a ). For example, with the holonomy corrections, 
we have \^g{y) + q{y) = z”/{k'^z) — kl{T]). Then, minimizing the error control function 


defined explicitly in (Zhu et al. 2014a), we find that in this case q{y) must be taken as 


g{y) = Once q{y) is determined, g{y) is in turn uniquely hxed. Then, the 

corresponding approximate analytical solution will depend on the number and nature (real 
or complex) of the roots of the equation g{y) =0 (Zhu et al. 2014a|[b|[c). In the quasi-de 

































Sitter background, it can be shown that g{y) currently has only one real root. In this case, 
the general expressions of the mode function, power spectra and spectral indexes up to the 


third-order approximations (in terms of A were given explicitly in (Zhu et ah 2014d). 


Applying them to the case with the holonomy corrections, we hnd (Zhu et al. 2015), 


~ A* 1 — 2 (l -|- Dp'j e*i — iA*e*2 + + 20* — 5 -|- 


'*1 


+ + ‘'*2+ 2 ^ 1 h ~ D * 5 ^ h ^*2 +(777“ + -^2 ) ^^*2^*3 


At 


TT" 1 

24 2^P 


a 


*2 


7 71' 


D* + — 7 + At -|- 2At j e*ie*2 — (^-D* -l- 6) 


A?(fc) = At 


1 + ~ 2 (^D* + 1 ) e*i — (4/1* + 6) 


TT 


2D*/ + 2D;+--5 + Al]et,+ 


-2D, 


*2 


TT' 


-Dp + — + 2At ) e*ie*2 


Us — 7 — 2e*i — e*2 + 4(5*/fe*i — 2e^i — (3 + 2D*) e*ie*2 ~ D*e*2f*35 


— —2e*i + 45*//e*i — 2e^;^ — 2 (/A* + 1) e*ie* 2 , — 16e*i(l +/l*e* 2 ), (8) 

wherein = P^/Pc -Cl, A^ = 181//*/(72e^7r^e*i), At = 181//*/(36e^7r^), D* = 67/181 —In 3, 
Zl* = 10/27 — In 3, At = ~ T’ At = sfiHs’ * denotes quantities evaluated at 

horizon crossing a(r 7 *)//(r 7 *) = ^/D{ri^k. e„ denote the slow-roll parameters, dehned as 
ei = —/////^, Cn+i = in/{Hen) {n > 1). Note that in the above expressions we have ignored 


terms at the orders higher than 0{e^^e^du)- To the hrst-order, it can be shown that our 


results are consistent with those presented in (Mielczarek 2014). 


In the case with the inverse-volume corrections, we have X^g{y) q{y) = k~‘^{z"/z — 
ojHv)), where ojl{g) is given by Eq.(|^, with Zs{g) = a(p[l + ^(ao — 2 ^ 0 )^^] and 
Zi,{g) = a(l — ao6pi/2), respectively. To minimize the errors, q{y) must be also chosen as 
in the last case, and then it can be shown that g{y) = 0 has only one real root, and as a 
result, the general expressions of the mode function, power spectra and spectral indexes 
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given in (Zhn et al. 2014d) are also applicable to this case, which yield (Zhu et ah 2015), 


1 — 2 (l + Dp) e*i — -D*e*2 + ^pi ( j f + Qq'' ^ 


^*(s) 




A?(fc) = A l-2(Zi; + l)e*i + ep; -i/J Q, 


•)*(*) 


ns — 1 — 2e*i — e *2 — 2e^]^ — (3 + 2Z1*) e*ie *2 — -D*e* 2 e *3 


+epi I -if* 




— ~2e*i — 2e^;^ — 2 (H* + 1) e*ie*2 + epi ( —H, 


K 


*{t) 
0 5 


r = 16e*i[l + iipe *2 - ep;(-if*)'^Q_i e. 




(9) 


Note that we parametrize Spiirf) = {api/ky{—ari)~'^y‘^ with epi = {api/ky, k = {—ar})~'^. 
In TABLE I, we list the values of the coefficients and for different 

values of a, as they represent the dominant contributions. The rest of the terms appearing 
in the above expressions are subdominant and will not be given here, but they are given 


explicitly in (Zhu et al. 2015). When a = 3, and vanish, so one has to consider 


contributions from and which are given by Qg^*^ = and = — |f'do- 

We emphasize that the modified power spectra and also spectral indices are now explicitly 
scale-dependent because of epi ~ 


Before considering the observational constraints, let us first note that in (Bojowald 


& Calcagni 2011 Bojowald et al. 2011a[b ) the observables n^, and r were calculated 
up to the hrst order of the slow-roll parameters. Comparing their results with 
ours, after writing all expressions in terms of the same set of parameters, say, 
ey[= M|,;(l/^/l/)V2], 7]v[= and we find that our 

results are different from theirs. A closer examination shows that this is mainly due to the 


following: (a) In (Bojowald & Calcagni 2011) the horizon crossing was taken as k = T-L. 


However, due to the quantum gravitational effects, the dispersion relation is modified to the 
form ([^, so the horizon crossing should be at Uk = ff. (b) In ([Bojowald fc Calcagni||20lT| 
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the mode function was first obtained at two limits, k ^ V, and /c -C and then matched 


together at the horizon crossing where k T-L. This may lead to huge errors (Joras & 


Marozzi| 2009 Ashoorioon et ah 2011), as neither nor is a good approximation 

of the mode function fik at the horizon crossing. The above arguments can be seen further 
by considering the exact solution of /i^. 


hfc(h)L=2 = 


Cl 


:WW 


zoi ly 




( 10 ) 


^/=r) ' ’' V 4^0^’ 2 ' 

for the case a = 2 , where hTlT( 6 i, 62 , denotes the WhittakerW function, oi = 1 — 
mepiK, 02 = 2j3oepiK, m{r]) is the coefficient of 5pi{r]) in the dehnition ^ + 

and V = 3/2 + 61+62/2 for the scalar perturbations, and v = 3/2 + 61 for the tensor. Matching 
it to the Bunch-Davies vacuum solution at k T-L, we hnd that Ci = e ^^I{y/2ka2 )• 
With the above mode function, the power spectra and spectral indexes can be calculated. 


and found to be the same as those given here, but are different from those of (Bojowald & 


Calcagni 

2011 ; 

Bojowald et ah 

2011 al|b 

). For more details, see ( 

Zhu et ah 

2015 


3. Detection of Quantum Gravitational Effects 


The contributions to the inflationary spectra and spectral indices from the holonomy 
corrections are introduced through the parameter which are of the order of 10 for 
typical values of the parameters ( Mielczarek||2014 ). Then, with the current and forthcoming 


observations (Abazajian et ah 2015), it is very difficult to detect such effects. 


On the other hand, for the inverse-volume corrections, let us consider the power-law 
potential E(0) = Xn4>"', for which we hnd that rjv = 2(n—l) 6 y/n, = 4(n —l)(n —2)6y/n^, 

where 6 y = Mp;n^/(20^). Thus, without the inverse-volume corrections {Spi = 0), we have 


ns{ev) and r = r( 6 v'), and up to the second-order of 6 ^, the relation (Creminelli et ah 
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2014), 


Tn{ns,r) = (n* - 1) + 


(2 + n)r 
8n 


(3n^ + 18n — A){ns — 1)^ 


= 0 , 


( 11 ) 


6(n + 2)2 

holds precisely. The results from Planck 2015 are rig = 0.968 ± 0.006 and ro .002 < 0.11(95% 


CL) (Komatsu et ah 2011 BICEP2/Keck and Planck Collaborations 2015 Planck 


Collaboration! 2015), which yields n < 1. In the forthcoming experiments, specially the 


Stage IV ones, the errors of the measurements on both Ug and r are < 10 ^ (Abazajian 


et ah 2015), which implies cT(r„,) < 10“^. On the other hand, when the inverse-volume 
corrections are taken into account {6pi % 0), we have Ug = ng{ey,epi) and r = r{ev,epi), 
and Eq. (0 is modified to. 


r„(n5,r) = J'(a) 


m 


( 12 ) 


where 6{k) = aoepiH'^ and J^{cr) ~ 0{1) (Zhu et al. 2015). Clearly, the right-hand side of 


the above equation represents the quantum gravitational effects from the inverse-volume 
corrections. If it is equal or greater than 0(10“^), these effects shall be well within the 
detection of the current or forthcoming experiments. It is interesting to note that the 


quantum gravitational effects are enhanced by an order Cy , which is absent in (Bojowald 


& Calcagni 2011) 


In the following, we run the Cosmological Monte Carlo (CosmoMC) code (Gong et 


al. 2008) with the Planck (Ade et al. 2013), BAO (Anderson et al. 2013), and Supernova 


Legacy Survey (Conley et al. 2011) data for the power-law potential with n = 1, which 


can be naturally realized in the axion monodromy inflation motived by string/M theory 


(Silverstein & Westphal 2008 McAllister et al. 2010). To compare our results with these 


acquired in (Bojowald et al. 2011b), we shall carry out our CMB likelihood analysis as 


closed to theirs as possible. In particular, we assume the flat cold dark matter model 
with the effective number of neutrinos A^e// = 3.046 and £x the total neutrino mass 
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Smj, = 0.06eV. We vary the seven parameters: (i) baryon density parameter, (ii) 

dark matter density parameter, (iii) the ratio of the sonnd horiozn to the angnlar 

diameter, 9, (iv) the reionization optical depth r, (v) 6{ko)/ev, (vi) ey, and (vii) Ag(/co)- 
We take the pivot wave number ko = 0.05 Mpc“^ used in Planck to constrain 6{ko) and 
ey. In Figj^ the constraints on d/ey and ey are given, respectively, for a = 1 and cr = 2. 
In particular, we hnd that 5{ko) < 6.8 x 10“^ (68% CL) for a = 1, and 5{ko) < 1.9 x 10“® 


(68% CL) for (j = 2, which are much tighter than those given in (Bojowald et ah 2011b). 


The upper bound for (5(fco) decreases dramatically as a increases (Bojowald et ah 2011b 


Zhu et al. 2015). However, for any given a, the best htting value of ey is about 10 which 


is rather robust in comparing with the case without the gravitational quantum effects (Ade 


et al. 2013). It is remarkable to note that, despite the tight constraints on S{ko), because of 
the enhancement of Eq.(12), such effects can be well within the range of the detection 


of the current and forthcoming cosmological experiments (Abazajian et al. 2015) for a ;< 1 


Note that small values of a are also favorable theoretically (Bojowald & Calcagni 2011). 


Conclusions 


Using the uniform asymptotic approximation method developed recently in (Zhu et 


al. 2014a|[b|[^ ), we have accurately computed the power spectra, spectral indices and the 
ratio r of the scalar and tensor perturbations of inflation in LQC to the second-order of 


the slow-roll parameters, after the corrections of the holonomy (Mielczarek 2008a, 2009 


Mielczarek et al. 2010 Grain et al. 2010 Cailleteau et al. 2012bpi Mielczarek 2014) and 
inverse-volume (Bojowald & Hossain 2007, 2008a|[b Bojowald et al. 2009, 2010 Amoros et 
al. 2014 Bojowald fc Calcagni 2011 Li Sz Zhu 2011 ; [Mielczarek et ah 2012 Bojowald et 


al. 2011b) are taken into account. The upper error bounds are < 0.15%, which is accurate 


enough for the current and forthcoming experiments (Abazajian et al. 


2015). Utilizing 
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the most accurate CMB, BAO and SN data currently available publicly (Ade et al. 2013 


Anderson et al. 2013 Conley et al. 2011), we have carried out the CMB likelihood analysis, 


and found constraints on {S{ko), ey), the tightest ones obtained so far in the literature. Even 
with such tight constraints, the quantum gravitational effects due to the inverse-volume 
corrections of LQC can be well within the range of the detection of the current and 


forthcoming cosmological experiments (Abazajian et al. 2015), provided that cr < 1. 


It should be noted that in our studies of the holonomy corrections, the effects of 
bouncing of the universe are insignihcant by implicitly assuming that inflation occurred 


long after the bouncing. This is the same as those considered in (Mielczarek 2008a, 2009 


Mielczarek et al. 1120101 IGrain et al. 1120101 ILi & Zhu||2011l IMielczarek et al. 1120121 ICailleteau 


et al. 2012bpi Mielczarek 2014). Thus, it is expected that quantum gravitational effects 
from these corrections are neglectible. However, when the whole process of the bouncing is 


properly taken into account, such effects may not be small at all (Grain & Barrau 2009 


Barrau & Grain 2014). It would be very interesting to reconsider the observational aspects 


of these effects, although cautions must be taken, as Eq.([^ was derived only for small 
potentials. Without this condition, there would be additional quantum corrections which 
are neither of holonomy nor of inverse-volume type, the condition. 


We thank Shinji Tsujikawa for valuable comments and suggestions. This work is 
supported in part by Ciencia Sem Fronteiras, No. A045/2013 CAPES, Brazil (A.W.); China 
NSF Grant, Nos. 11375153 (A.W.), 11047008 (T.Z.), 11105120(T.Z.) and 11205133 (T.Z.); 
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Table 1: Values of Coefficients and for different values of a. 
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Fig. 1.— Two-dimensional marginalized distribution for the parameters S/ey and ey at the 
pivot ko = 0.002Mpc“^ for the power-law potential with n = 1. (a) The left panel is for 
a = 1, and (b) the right panel is for <7 = 2. The internal and external lines correspond to 
the conhdence levels of 68% and 95%, respectively. 





